tic  UTi-cwm  ytesLEK  a*©  tte 
K5SIiMXnm>3nEFEL  TCASSFGWSATIOS; 


C-E-O.  Ciccsflia 

The  Sslftrsltf  oi  Texas  a  Asst in 
Assets,  Texas 


AMSL  1.018 
August  1979 


Applied  Mechanics  Research  Laboratory 
The  University  of  Texas  at  Austin 


Austin,  Texas 


This  report  was  prepared  under 


T&Acr 

-£*■♦**  Bo.  *000  U-$7A-0126-0013 

for  the 

Office  of  Havel  Research 

by  the 


Applied  Mechanics  Research  Laboratory 

The  University  of  Texas  at  Austin 
Austin,  Texas 


under  the  direction  of 
G.E.O.  Giacaglia 

Visiting  Professor,  Aerospace  Engineering 
and  Engineering  Mechanics 


THE  LEVI-CIVITA  PROBLEM  AND  THE 
KUSTAANHEIMO-STIEFEL  TRANSFORMATION 


by 


G.  E.  0.  Giacaglia* 

The  University  of  Texas  at  Austin 


*0n  leave  from  the  University  of  SSo  Paulo  and  the  University 
of  Campinas,  SSo  Paulo,  Brasil. 


THE  LEVI-CIVITA  PROBLEM  AND  THE 
KUSTAANHEIMO-STIEFEL  TRANSFORMATION 

G.  E.  0.  Giacaglia* 

The  University  of  Texas  at 
Austin 


Abstract 

- >  By  imposing  regularity  of  the  new  accelerations  and 

»  , 

velocities  we  introduce^  a  time  and  a  coordinate  transforma¬ 
tion  which  generate  as  particular  cases  the  Levi-Civita  and 
the  KS-transformations . 


Introduction 


Consider  the  differential  equation 


,3w.t 


(1) 


where  x  is  an  n-column  in  R  ,  W  *  W(x)  is  a  scalar  function 

and  3W/3x  its  row  gradient  with  respect  to  x»  The  superscript 

T  indicates  the  transpose  of  a  vector  or  matrix.  Let  x  be 

n 


“  2  % 

the  absolute  value  of  x,  x  *  (}  x.)  . 

J-l  3 


We  consider  the  case 


W  -  1 / x  +  R(x)  +  E 


(2) 


and  assume  that  R(x)  is  in  an  open  region  ft  of  Rn  containing 
the  origin  x  =  0.  It  follows  that  for  every  x  in  ft  except 
x  *  and  for  every  finite  x  there  exists  and  is  unique  a 


*0n  leave  from  the  University  of  Sao  Paulo  and  the  University 
of  Campinas,  Sao  Paulo,  Brasil. 
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solution  of  Equation  (1)  which  can  be  continued  in  ft  except 
through  the  origin.  We  wish  to  study  the  solutions  which 
approach  the  singularity  x  -  0.  Assuming  E  to  be  an  arbi¬ 
trary  constant,  the  energy  integral  can  be  written 

xTx  *  2W(x)  (3) 

Let  cs  consider  the  tangential  transformation 

dx  =  2M(u)  du  (4) 

where  u  is  an  n-column  in  Rn  and  M  an  n  x  n  matrix.  Necessary 
properties  are  that  in  an  open  region  fty  containing  the  origin 
u  «  0.  and,  except  at  a  finite  number  of  points,  M  ^  exists 
and  also 


(5) 


for  k,£  =  l,2,...,n  and  s  >  £,  which  are  the  integrability 
conditions  of  Equation  (4). 

We  also  introduce  a  new  independent  variable  t  by 

dT  e  a(u)  dt  (6) 

2 

where  a  (u)  is  a  scalar  function,  at  least  C  in  flu  except  at 
a  finite  number  of  points.  With  the  use  of  Equations  (4)  and 
(5)  into  Equations  (1)  and  (3),  we  obtain  (Giacaglia,  1968) 
the  following  relations 


u”  -  -M-1Mu*  -  <o*/a) u*  +  (l/4o2)  M^OT1)*  (3V/3u)T  (7) 

and 

u’TMTMu’  -  W/2o2  (8) 

where  prices  indicate  derivatives  with  respect  to  t. 

Regularization 

We  require  that  at  x  *  0,  |u£|  <  A,  A  finite.  Ve  also 
restrict  ourselves  to  orthogonal  transformations ,  that  is, 
such  that 

g(u)<J2(u)MTM  -  I  (9) 

o 

where  I  is  the  identity  n  x  n  matrix  and  g(u)  is  a  scalar  C“ 
function  in  fiu  except  at  a  finite  number  of  points.  Equation  (8) 
becomes 

u’2  -  |gW  (10) 

while  (7)  can  be  written  as 

u"  =  F(u,u')  +  i  Cf^(gW)]T  (11) 

where 

_F(u,u'  )  =  -go2MTM'u'  -(a'/o) u'  -  (l/2g)  (3g/3u)Tu'  xu* 

(12) 

It  follows  from  (10)  that  if  |  |  <A  as  x  0,  then  g  must  be 

infinitesimal  of  order  x,  that  is,  in  the  neighborhood  of 
x  =  0  we  may  write  g  *  xf (x)  where  f (£)  is  finite  and  not  zero. 
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For  the  transformation  to  be  canonical  it  is  now  sufficient 
to  have  *  «  0  whatever  u* .  The  p-th  component  of  F  is 


F  *  -go  T  n.  C?M.  /3b  }u'u*  -  —  Toe/au  )u  u* 
P  °  .  L  Ip  is  CC*  O  a  a  p 

r  j<t,q  q  ’ 


-(l/2g)  (3g/3u  )  Jo*2 
F  q  q 


For  ?  to  he  identically  zero,  the  coefficients  of 
? 

(u'u*)  must  vanish  whatever  (a, 6)  are,  so  that  we  obtain 
a  p 


3M. 


2  2  J  riB  +  l(|o_4  +  |*_6  )ti  at<  , 

8  "  jp  »o  a  3u  gp  3n.  op  g  cu  ag 

3  a  a  S  * 


(13/ 


for  all  values  of  p,a,g.  If  one  multiplies  Equation  (IS) 
by  (M_1)  .  and  nakee  use  of  Equation  (9)  the  result  is 

pi 


a*. 


■Aft.  4.  A —  iSLOL —  u  4  -  u  ) 

3u  2 j  *3u  ig  3u  nicr 

a  a  g 


+  h  «„«  IHi  -  0 


2g  og  £  ip  3up 


(14) 


for  all  (i,a,g),  the  equation  being  symmetric  wit-h  respect 

2 

to  the  pair  (a,g).  This  is  a  system,  therefore,  of  n  (a+l)/2 
partial  differential  equations  which  define  K,  g,o  such  that 
the  vector  F  is  identically  zero. 

Getting  a  new  system  independent  of  the  velocity  u1 
and  performing  an  orthogonal  transformation  will  reduce  the 
equations  to  the  canonical  form 


5 


I  r£_ 
4  L3u 


(«W)]T 


(15) 


and  this  can  be  Versed  to  be  Levi-Cirita* a  problea.  In  pria- 
clple,  the  probier  is  now  reduced  Co  che  choice  of  Che  scalar 
function  g(u). 

Conneccion  between  cine  and  coordinates  transformations 

A  si»j»2.e  relation  between  g  and  a  can  be  obtained  froa 
Equation  (13).  In  fact,  considering  the  case  0  •■=*  p  #  a,  one 
finds 


,  2  v  3M.  .  . 

2go  1  M  _ IE  +  ±  l2_ 

3  jp  3«a  a  au^ 


or 


1  JL<L_  ,  1_  l  m2 
goJ  3u  3u  jp 


or,  considering  Equation  (9), 


1  3  o 
go*  3u 


3  .  1  . 

TT  W' 


so  Chat 


3u 


(in  go) 


for  all  values  of  a.  Thus,  by  proper  normalisation 


go 


(16) 


independently  of  the  form  of  matrix 


out  its  orthogonality. 
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This  result:  indicates  the  fact  that  if  g  Is  infinitesimal 
with  x  as  resulted  from  Equation  (10) ,  then  1  /o  defining  the 
time  transformation  must  have  exactly  the  same  character  as 
x  -*■  0.  It  also  follows  from  Equation  (9)  that  jdet  Mj  *  |g|n^, 
so  that  if  and  only  if  n  *  2  the  jdet  Hf  is  infinitesimal 
of  the  same  order  of  g  as  x  0. 


The  choice  of  |(u(x))  and  the  connection  between  vi  and  x 
Use  of  (16)  into  (14)  gives 


16  _  1_ 

3u  2g 

o  * 


(M  -  +  M  -  S  .  [  H.  )  =  0  (17) 

iS3u_  io3u_  ag  L  xp3u 
"  P  P 


The  sinplest  choice  of  g,  satisfying  the  above  derived 
regularization  condition  is 


g  =  Kx 


where  K  is  a  constant,  not  zero. 


It  follows  that 


I2-  — -  *  2K l  —  M 

3u  *r  3x.  3u  y  x  1 

p  i  2.  p  £ 


H7  83u  J  *2  M£p 

p  £ 


Multiplying  this  last  by  (M  )  and  taking  Equation  (9)  into 

P* 

account,  summation  of  p  gives 


! 

i 


7 


fa_ 

Su 

P 


x 

a 


(18) 


Differentiating  again  vice  respect  to  and  multiplying  by 
X  ,  ve  obtain,  after  summation  over  a.  and  together  with  the 
Sntegc ability  condition. 


(19) 


This  last  set  of  n  partial  differential  equations 

defines  the  function  g(u>  *  Kx,  explicitely  in  teras  of 

u, ,  u_,...,u  .  Evidently  this  is  cot  enough  to  define  the 
1  i  n 

transformation  froa  x  to  ti,  or,  better  say,  the  matrix  M. 

Up  to  this  aoaent  one  can  conclude  that  for  an  ortho- 

T 

gonal  and  ictegrable  matrix  ft(u),  such  that  M  M  =  gl,  if 
g  *  Kx  and  the  new  equations  are  independent  of  the  velocity 
u‘ ,  then  the  possible  transformations  M  should  be  such  that 
g(u^)  satisfies  Equation  (19).  Some  of  the  results  of  these 
sections  were  presented  in  a  previous  work  (Ciacaglia,  1968) 
related  Co  the  n-body  problem. 

Two-dimensional  case.  Levi-Civita's  Transformation. 

The  orthogonality  condition 
MTM  =  gl 

can  be  written,  in  general,  as 


8 


*11 

*  /S 

cos 

*12 

-  /£ 

sin 

*21 

-  ^ 

sic 

*22 

» 

sin 

(20) 


where  $  *  $(u)  has  to  be  determined. 

By  considering  the  integrability  conditions 

3*ll^3u2  “  3*12^3ul  *  3M21^3u2  *  3*22^3uI  *  we  obtain 


14.  ,  1_  Is.  M  ,  1& 

Suj  ~  2g  3u2  *  3u?  2g  3U]l 


2  2 

or,  since  one  muct  have  3  $/3u^3u2  =  3  $/3u23u^,  the  result  is 


2g  & 


i  L  3ui 


>2  +  (|*  >2 

3u2 


(21) 


But,  for  n  =  2,  Equation  (19)  gives 


2g  >uf  + 


rM  }2  +  (l£.  }2 

u.  3u0 

X  4. 


-  8K2g 


so  tnat,  necessarily. 


M  -  2K‘ 
2 “1 


and  similarly 
2 


f-f  =  2iT  . 

3u| 


By  making  use  of  Equation  (19)  we  finally  obtain 


2  2 

..1.2.  2*  .  .  .  a  +  8 

g  *  K  (Uj  +  «2)  +  auj  +  8u2  +  — - 


with  a , B  arbitrary  constants.  It  is  worth  noting  that  the 
following  solutions  of  Equation  (21) 

g  *  -|  K2  (Uj  +  u2^  +  2aK^~  ^ul  ±  u2^  +  “2 

do  not  lead  to  orthogonality  and  integrability  of  M  so  they 
must  be  discarded.  There  is  no  other  possible  solution. 

For  the  angle  4  we  obtain 


M 

3u, 


<2K  u2  +  8) 


&  *  +  H  <2k2“i  +  “> 


or,  integrating. 


=  arctan 


2K  u2  +  g 


•uL  +  a  '  Y 


where  y  is  an  arbitrary  constant.  We  now  obtain 


/g  cos  4>  =  [(2K2Uj  +  a)  cos  y  +  (2K2u2  +  3)  siny]/2K 
/g  sin  <j>  =  [-(2K2u1  +  a)  sin  y  +  (2K2u2  +  3)  cosy]  /2K 


so  that  the  transformation  matrix  M  can  be  written 


M  =  R(y)  K 


U1  "u2 
u2  U1 


a  -8 
+  8  a 


where 
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R(y>  -  f  TY  slnY| 

l~sinY  cos yj 

la  complex  fora,  the  transformation  can  be  written 

dx*  *  2eilf(Ku*  +  a*)du*  (24) 

where 

x*  *  Xj  +  ix^  ,  a*  *  Uj  +  iUj  ,  a*  *  a  +  iB  » 
Equation  (24)  gives  the  most  general  form  of  Levi-Civita’ s 
transformation  which  corresponds  to  the  particular  case 
I  ■  0,  K  »  1,  o  ■  B  *  0.  In  the  same  notation  of  Equation 
(24)  the  function  g  is 

-  g  «  (Kut  +  a)2  +  (Ku2  +  8)2  (25) 

with  the  corresponding  Levi-Civita  case 

2  2 

8  =  U1  +  u2  ' 

We  also  remember  that  the  transformation  is 


dx  =  2Hd 


so  that  one  finds 


2  21 

’  ' 

““ 

u,  -  u„ 

a  -  6 

x  =  R(y) 

K 

1  z 

2ulu2 

+  2 

6  a 

(26) 


or,  in  complex  notation 

x*  =  e^Y(Ku*2  +  a*u*)  +  b*  , 
where  l)  is  an  arbitrary  vector.  It  can  also  be  written  as 


11 


obtain 


12 


3f  33h  .  33f  3h  Y  32f  3zh 

3x  Tp  TP  3y  3x3  u  3y3u 

a  a  3  a 


Admitting  that  (3f /3x) • (3f /3y)  t  0,  division  by  this 
quanti^v  leads  to 


3  3f / 3x3  33h/3v3 

3£/3x  3h/3y 


K(x,u  )+  L(y,u  )  = 

Cl  vt 


3  2f  32h  1 

\  !\ 

3f 

3h 

3x3u  3y3u 

a  a 

) 

l/i 

3x 

3y 

A 

Again,  successive  differentiation  with  respect  to  x 
and  y,  gives 


f  3h  33h  3  2h  32h] 

3 f  3  3f  32f  3 2 f 1 

3y  3y^3u  3y3u  Tp 

3x  3x2au  3x3u  TP 

a  a 

a  a 

and  this  relation  has  to  be  valid  for  all  pairs  (k,£),  that 

is,  whatever  u  is. 

a 

The  only  way  to  satisfy  both  the  last  equation  and 
Equation  (19)  is  then  given  by  the  2n  equations 


3 h 

3  3h 

9  2h 

3  2h 

9y 

3y^3u 

a 

3y3u 

a 

TP 

31 

3  3  f 

3  2  f 

.  2  . 

3  r 

3  x 

3x23u 

a 

3  x3u 

a 

3  x  2 

more 

specifically, 

3  2  f 

TP 

2A  = 

const 

3  2h 

TP  = 

2P  = 

const. 
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for  every  pair  ( k,l ).  It  follows  from 
f  =  Ax2  +  B(uo,y)x  -f  C(ua,y) 
h  =  Py2  +  Q(ua,x)y  +  R(ua,x) 

that  the  most  general  allowed  form  for  g,  under  the  above 
specified  conditions,  is  given  by 

n  9 

g  *  I  <Ku  +  a  r  (30) 

j=l  2  2 

where  K  and  the  a^  are  constants,  and  there  remains  to  be 
verified  whether  this  form  leads  to  integrable  relations  of 
the  matrix  M.  More  details  on  these  results  were  given  by 
Giacaglia  (1966). 

The  definition  of  the  transformation  M 

Consider  again  Equation  (18)  and  differentiate  it  with 
respect  to 


Let  us  define 


14 


to  *  Ku  +  a 
P  P  P 


(p  *  1,2,. • • , n) 


(32) 


so  that  Equation  (31)  becomes 


9M 


mi 


p  9  to 
P  P 


M 


mi 


and  therefore  is  an  homogeneous  function  of  the  first 

degree  of  the  to  ,  that  is,  a  linear  function  of  the  u  . 

P  p 


Let 


Mk*  *  l  a 


kius 


(33) 


where  the  matrix  M  must  be  orthogonal  and  lead  to  integrable 
relations . 

Consider  first  the  orthogonality 


l  "kAj  -  *  sii 


where,  according  to  (30)  and  (32) 


- 1  • 


(34) 


We  obtain 


y  l  (  la,  . ai’.)to  to  *■=  6..  T 

s  r  kkik3  sr  ij  J 


to 


or 


s  r  r  s  0 .  . 

akiakj  +  akiakj  "  26ijdrs 


Consider  the  n  matrices 


(35) 
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so  that  Equation  (35)  can  be  written 


T  T 

aa  +  a  a  *  21  6 
s  r  r  s  rs 

with  the  integrabllity  conditions  giving 


(36) 


P  ! 

Ctf.  “  O. 

k£  kp 


(37) 


The  solution  for  u  =  2  is  given  by 


1  0 
0  1 


f°  -l] 

a2  “  [l  Oj 


giving  back  the  general  Levi-Civita  transformation 


M 


U), 


’V 


(0 


1] 


For  n  =■  3,  there  is  no  solution,  while  the  problem  can  partially 
be  solved  for  n  =  4  or  n  =  8  (Hurwitz,  1933). 

In  order  to  study  the  problem  for  n  =  4  we  here  make 
use  of  quaternions  as  opposed  to  the  above  spinor  solution 
of  the  case  n  =  2. 


The  four-dimensional  case 

We  follow  the  results  in  Hurwitz  (1919)  and  define  the 
imaginary  units  i,j,k,  satisfying  the  multiplication  rules: 
i2  =  j2  =  k2  =  -1,  i-j  -  -j*i  =  k,  etc. 

The  rotation  matrix  M  is  then  generated  by 

dx^  +  idx2  +  jdx3  +  kdx^  =  2(f1  +  if2  +  jf3  +  kf^)* 

•  (du^  +  idu2  +  jdu3  +  kdu^) 
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which  we  write  dx* 
represented  by 

X 

f. 


f  3 


2f*  •  du*.  We  easily  find  that  M  is 


-f. 


-f . 


-f.  -f , 


-fi 


-f. 


(38) 


where,  as  was  seen,  the  f^  must  be  linear  functions  of  the 
Uj  or  homogeneous  of  first  degree  of  the  .  The  above  matrix 
is  orthogonal  whatever  the  f^  are  (not  all  zero).  The  matrices 
otj  (j-1,2,3,4)  are  obtained  immediately  by  setting  one  of  the  f ^ 
equal  to  unity  and  the  others  equal  to  zero,  in  succession. 

By  imposing  integrability  conditions  on  the  first  three 
rows  of  M  we  obtain  the  conditions 


fl  “  *^ui» 

U2.Uj) 

^2  =  ^2^U1 

»u2»u4> 

f3  =  f3(ul 

*u3  »u4> 

f4  '  £4Cu2 

,U3 ,U4> 

where 

<{>  is 

linear 

while  f2, 

f  3  ’  f  4 

must  satd 

Lsfy  the 

'  3f2 

/dUj  = 

—  3 4> /  a u2  , 

3f2/9u2 

=  3<j>/3u  ^ 

,  3f2/3u 

(3g)  - 

af3 

/  SUj,  = 

— a  <j>  /  a  u  ^  » 

3f3/3u3 

=  a  /  a  u  j 

,  3f3/3u 

3  f  4 

/9u2  = 

— a  4> /  a » 

3f .  /3u„ 

4  o 

=  3  <j>  /  a  w  2 

,  3f4/3u 

If  integrability  of  the  fourth  row  is  imposed,  there 
results  <j>  =  0,  so  that  complete  integrability  is  impossible. 
In  other  words,  whatever  d>  is  (^0),  the  differential  form 


f^duj  -  fjdu^  +  f?du^  +  f jda^  is  not  integlable.  This  prob- 
iei,  of  course,  does  not  represent  an  obstacle  but  only  means 
chat  to  the  holocomic  constraint  i,  »  0  fo  the  physical  space 
there  corresponds  the  con  holononie  constraint  (£*  *  du*)^  ** 
in  the  u-spacc.  Tbis  problem  has  been  discussed  in  several 


works  for  the  particular  transformation  (KS-Traes formation) 
introduced  ly  Kustaanheimo  and  Stiefel  (1565}. 


We  cow  obtain  the  general  solution  of  the  problem  by 


setting 


fj  *  ®i“i  +  ®2u2  +  S3U3 


so  that,  applying  Equation  (39),  we  easily  find 


f2  -  -62»,  +  8j»2  -  6jU4 


f3  *  -e3wl  +  61m3  + 


where 


f4  '  -s3“2  +  V3  '  6l“4 


u .  =  Ku .  +  a . 
J  J  3 


The  transforaation  M  depends  therefore  on  seven  para- 
aeters,  while  K  is  just  a  scaling  factor  which  we  can  set 
equal  to  unity  without  loosing  anything*  Sot  all  B's  can  be 


In  components  form,  the  transformation  can  now  be 


written 


18 

EdXj  *  3jd(asj  -  acj  -  +  + 

-r  -  ae^/  + 

+  S,d(UjUj  +  a^M^) 

Kdx2  «  Sjdtej*^  +  u3«4)  + 

+  89d(-wj  +  Mj  -  «j  + 

+  p^d(-ajOS^  + 

Edx,  *  SjdfMj^  -  «2e4^  + 

+  E9d(Mj«4  +  a *2“^  + 

A  e  .,  2  2  2  2. 

+  S3d(-Uj  -  «2  +  Uj  +  «i) 

Kd*^  *  Sj(-e*4daij  -  tt^dca,  4-  »,d«3  +  Bjde^)  + 

+  B2(«,dUj  -  B^iu^  -  BjdtSj  +  OjdM^)  + 

+  g,(-»_do,  +  as, da.  -  u.du.  +  e>,dfe,)  (42) 

321  12  43  34 
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«4.) 


and  it  is  clear  that  dx^  is  not  an  exact  differential.  Setting 

K  =  1,  a.  *  0  (j=l, 2, 3, 4),  8,  =  S,  =  0,  we  obtain  essentially 
j  z  a 

the  KS-Transforaation. 

Other  than  transf orcation  (42)  one  could  consider  rota¬ 
tions,  translations  and  inversions  which  are  analytic  and  can 
be  used  to  Bake  the  transformation  satisfy  other  conditions 
than  the  ones  required  in  this  work.  For  example,  they  were 
used  by  Waldvogel  (1967)  to  generalize  Birkhoff's  transforma¬ 
tion  for  the  3-dimensional  restricted  problem  of  three  bodies. 
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In  that  work,  Che  anchor  gives  as  excellent  discussion  of 
the  general  problem. 

Generalised  ferns  of  generators  for  regularization 
were  given  by  Szebebely  £19i7a),  Giacaglia  (1967a),  Giacaglia 
(11966),  Giacaglia  (!967b),  and  they  are  totally  discussed  by 
Szebehclv  (x96?b). 

Tbe  Jesuits  given  by  (42)  show  that  it  is  not  possible 
to  fully  generalize  Equation  (24)  into  a  four-dimensional 
complex  analogy  (quaternion  formalism)  which  would  look  like 

dx*  =  2Q*(S)  •  (Ku*  +  a*)  -  d>u* 

with  Q*(5)  a  rotation  quaternion.  It  is  easy  to  verify 
that  the  quaternion  product  u*  •  da*  is  not  intcgrable.  The 
KS-Transforaation  corresponds  to 

dx*  *  2(u,  +  iu2  +  ju^  -  ku^MdUj  +  idu2  +  jdu^  +  kdc^> 
=  2u*  -  du*  ,  (43) 

where  u*  =  u1  +  iu2  +  ju^  -  ku^  , 
while  (42)  can  be  written  as 

Kdx*  =  2(3ju£  •  du*  -  &'£a%  *  ”  ®3UZ  *  (44) 

with  the  same  notation  of  Equation  (43),  that  is, 

u*  =  Uj  +  iu2  -  ju3  +  ku^ 

u*  =  Ui  -  it  2  -r  ju3  +  ku4 


All  combinations  u£  *  du*  ,  u|  •  du*  ,  ui|  •  du*  cor: aspond 
to  integrabic  relations  leading  to  orthogonal  transformations. 
In  any  event  the  tine  transformation  is  given  by  Equation  (34), 
that  is. 


where  as  »  Ku^  +  a.  ,  an  exact  generalization  of  Levi- 
Civita's  general  correspondent  Equation  (25). 

It  is  important  to  notice  that,  in  quaternion  notation, 
the  integrated  (3  ■*  4)  ES-Transformation  can  be  easily  written 
as 


or 


xt  +  ix2  +  jx3  -rkd  *  (Uj  +  iu2  +  ju3  +  ku^)  * 

*  (u3  +  iu2  +  ju3  -  k»4) 


X*  =  u*  •  u£ 


In  the  same  sannet  one  can  write  the  integrated  form 
of  (44)  as 


Ks*  *  6,u*  *  wf  -  S«u*  *  u*  -  g_u*  •  u* 

1  4  2  3  3  2 

the  k-th  component  of  this  being  identically  zero. 

Quaternion  notation  has  been  neglected  for  a  long 
time  but  just  recently^  Arenstorf  (1969/  has  shown  its  use¬ 
fulness  in  relation  to  the  transformation  introduced  by 
Waldvogel  (1967). 
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Conclusions 

The  primary  reason  for  this  paper  was  to  obtain  KS- 
Transformation  via  the  direct  solution  of  the  original 
problem.  As  a  result  we  obtained  a  system  of  Partial 
Differential  Equations  from  which  arose  the  equation  con¬ 
necting  the  time  and  coordinates  transformations.  We  then 
obtained  a  general  form  for  the  time  transformation  and  with 
this  we  concluded  that  the  coordinates  transformation  matrix 
had  to  be  linear  in  the  new  variables,  or  homogeneous  of 
first  degree  in  a  veil  defined  linear  combination  of  these. 
These  combined  conditions  led  to  a  transformation  which 
contains  the  KS  as  a  special  case.  It  is  important  to  note 
that  orthogonality  of  the  Jacobian  matrix  for  the  transforma¬ 
tion  and  absence  of  velocity  dependent  terms  was  imposed  as 
a  starting  point. 
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